Escaping from nonhyperbolic chaotic attractors 
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The noise-induced escape process from a nonhyperbolic chaotic attractor is of physical and 
fundamental importance. We address this problem by uncovering the general mechanism of 
escape in the relevant low noise limit using the Hamiltonian theory of large fluctuations and by 
establishing the crucial role of the primary homoclinic tangency closest to the basin boundary 
in the dynamical process. In order to demonstrate that, we provide an unambiguous solution of 
the variational equations from the Hamiltonian theory. Our results are substantiated with the 
help of physical and dynamical paradigms, such as the Henon and the Ikeda maps. It is fur- 
ther pointed out that our findings should be valid for driven flow systems and for experimental data. 
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Many nonequilibrium systems in physics, chemistry, 
biology or technology exhibit, as a crucial feature, noise- 
induced escape from a metastable state. Examples range 
from Josephson junctions 0], switching in lasers 0], Pen- 
ning traps 3], over chemical reactions '3| and protein 
folding ^] to electronic circuits 6, 7]. For small noise in- 
tensities, a WKB-like extension of Kramers' equilibrium 
theory ^ has been developed to treat the realm of 
nonequilibrium systems ,10. llj. This approach, mak- 
ing use of an auxiliary Hamiltonian system, identifies the 
so-called most probable exit path (MPEP), whose proba- 
bility of occurrence is exponentially larger than that of 
all the other paths. The MPEP can be observed through 
its 'prehistory probability distribution' and it was carried 
out numerically (13i | and experimentially 

The variational equations for the calculation of the 
MPEP are well known both for continous 0, ^| and 
discrete systems 0, 0, 0, ^3 ' Y^^ it is in general not 
clear how to solve them. Methods have been proposed 
only for the special case of escaping from periodic states 
(fixed points, periodic orbits or limit cycles) H^ . 
Thus the challenging question of noise-induced escape 
from a chaotic attractor (CA), although tried before, re- 
mained to be answered. The reason is that all previous 
work dealing with the escape from a CA relies on Monte 
Carlo simulations of the escape path. By making use 
of the 'prehistory probability distribution', an unstable 
periodic orbit embedded in the CA was identified and 
taken as initial condition for the calculation of the MPEP 

In this Letter we show that, in fact, the initial condi- 
tion for the MPEP, and thus the path itself, is uniquely 
determined by the primary homoclinic tangency (PHT) 
closest to the basin boundary, as well as its preimages 
and images. Our solution of this long-standing problem 
does not only uncover the general mechanism of noise- 
induced escape from a nonhyperbolic CA, but it reme- 



dies, in addition, shortcomings in the previous method. 
Herewith, we establish the following. First, no unsta- 
ble periodic orbits have to be determined and no Monte 
Carlo simulations have to be run to identify these orbits 
in the escape path. Second, the arbitrariness of the selec- 
tion of one particular periodic orbit, whose coincidence 
with the MPEP is not given from first principles, can be 
dispensed with. As a matter of fact, the PHT is, in con- 
tradistinction, a completetly deterministic quantity and 
can be both easily calculated numerically [23j and ex- 
tracted from experiments (241] . Third, complementary to 
the work of Ref. 22] in which the importance of the ho- 
moclinic structure of the fractal basin boundary to the 
MPEP is stressed, the present work closes an important 
gap in pointing out the essential role of the homoclinic 
structure of the chaotic attractor. 

Virtually all CAs occurring in nature or serving as 
prototype models in nonlinear dynamics are nonhyper- 
bolic. Nonhyperbolicity means that the stable and un- 
stable manifolds of the system are tangent in the phase 
space. If both manifolds belong to the same periodic or- 
bit, these tangencies are homoclinic. They are called pri- 
mary, if the, generally quadratic, curvature of the man- 
ifolds in the vicinity of the tangent points attain a min- 
imum. Bounded noise on nonhyperbolic CAs causes at- 
tractor deformations |23ll24| . which are most pronounced 
at the forward and backward iterations of the PHT clos- 
est to the basin boundary. As a consequence, PHTs can 
be better defined as being those tangencies which ex- 
hibit an amplification of a pertubation under forward and 
backward dynamics. 

In the following we consider discrete systems, but our 
results are not limited to them and apply equally well 
to any flow that can be reduced to a Poincare map. For 
the nonequilibrium case, in analogy to Kramers' law, the 
mean first exit time is given through the least action S 
[ilESEi by (r) - exp [■§] , where D is the variance of 
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the additive Gaussian white noise. To be specific, for a 
d-dimensional map Xn_|_i = f(xn) + the action of the 
escape trajectory to be minimized has the form 



1 ^ 



ffx. 



1 ^ 

9 ^ n , 



(1) 



where the are the noise vector terms. The boundary 
conditions are such that Xi is a point of the attractor and 
xn is on the basin boundary, from where no fluctuations 
are needed to pass to another stable state. The MPEP, 
which minimizes this action, can be calculated through 
the Lagrangian ,12. .16.] 

N N 

i= 2E^n^n + E-^n(Xn+l-f(Xn)-en) (2) 

n— 1 n— 1 

to yield upon variation of ^n,An, and x„ the area- 
preserving equations 

x„+i = f(x„) + A„ (3) 
A„+i = {(Df(x„+i))T}-^ An, (4) 

where Df is the Jacobian matrix of f . The Lagrange 
multipliers An replace the noise terms ^n- The opti- 
mal solution of Eqs. 13I4|I yields the least action 5* = 
lS^=i-^n -^n, the corresponding MPEP (given by the 
Xn) and the optimal force (An). 

The solution of these equations is intricate though, 
caused by wild, fractal fluctuations of the energy land- 
scape having many local minima and maxima 25]. A 
way to deal with these difflculties, when escaping from a 
periodic solution, is to employ a refined shooting method 
|l8l Hgf . It consists of a parameterization of the initial 
conditions Xi in Eq. Q on a small circle centered in 
one of the components of the periodic orbit Xpo , and for 
Ai pointing to the unstable manifold of Xpo- Equiva- 
lently, one can take as an initial condition one point of 
the periodic orbit itself xi = xpo and Xi G r x cj), where 
(j) G [0, 2 tt] and r is varied within r e [rmin, L Tmin] ■ Here 
L is the largest eigenvalue at x po and Tmin arbitrary, but 
small. The upper limit of r guarantees that every solu- 
tion is only considered once. This provides an efficient 
numerical method for the calculation of the MPEP. 

Now, we show that the same procedure can be applied 
to calculate the escape from a chaotic attractor. To do 
this, one has to select the initial condition xi as a preim- 
age of the PHT closest to the basin boundary, since there 
noise causes the largest deviation This means that 
the energy to leave the CA is the lowest. The param- 
eterization for Ai remains unchanged, since every point 
on the CA has a well defined largest eigenvalue. We first 
demonstrate this for one of the funamental dynamical 



paradigms, the Henon map |2g 

yn-\-l • 



(5) 



We choose the parameters a = 1.3 and b = 0.29 for which 
a CA and another attracting state at infinity coexist. In 
Fig. ^the CA is plotted, as well as its basin of attraction, 
and the saddle point on the boundary (square). We also 
include the PHT closest to the boundary together with 
its 10- fold preimages (circles) and the MPEP (crosses). 
The latter is found by iterating Eq. at the 10-fold 
preimage of the PHT and by looking for the absolute 
minimum when changing Ai as described above. As it 
can be seen, the MPEP moves initially very close to the 
deterministic dynamics yet it deviates increasingly with 
every iteration (as expected, since otherwise there were 
no escape). At the PHT the MPEP already differs con- 




-2.5 



-0.5 



1.5 



FIG. 1: (color online). CA for the Henon map with a = 1.3 
and b = 0.29. The hatched region markes the basin at infinity. 
The PHT at {x,y) = (1.653, -0.099) is shown, together with 
its 10-fold preimages (numbered -1,..,-10) as circles, the 3-fold 
images of the PHT (numbered 1,2,3) as diamonds, while the 
MPEP is depicted with crosses, connected through a dashed 
line to guide the eyes. The saddle point on the boundary is 
marked with a square. 



siderably from it. Only there, it is for the first time 
located outside the CA and advances then, following the 
elongated images of the PHT (depicted as diamonds), to 
the basin boundary. From there it approaches the saddle 
point located on the basin boundary. We emphasize that 
it is not presupposed that the MPEP leaves the CA in 
the vicinity of the PHT. We simply use a preimage of the 
PHT as the initial condition Xi and look for the optimal 
solution of Eqs. l|2E|l. The fact that the MPEP is indeed 
following very closely the structure of the PHT (preim- 
ages and images) is thus a confirmation of our claim. 

To shed further light on this general mechanism, the 
stepwise action S'„ = ^Aj An is plotted in Fig. ^ There 
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FIG. 2: The stepwise action Sn ~ fA^ An for the Henon map. 
At each time step, the value of the graph corresponds to the 
optimal force An, as a result of the optimal solution of Eqs. 
(I3I4II . The inset shows the data in semi-log scale. 



is only one large peak (step 10) where the MPEP moves 
away from the CA close to the PHT. Afterwards, the re- 
quired energy drops off exponentially, since the MPEP 
follows a relaxational path. It is important to note that 
this scenario differs from the one of the attractor defor- 
mation. In the latter, the trajectory deviates from the 
CA starting at the PHT, while in our case the MPEP 
departs already at the preimages (cf. Fig. |21 inset) and 
is, when passing the PHT, a finite distance away. This 
is the reason why we have to start iterating Eq. ^ at 
some preimage of the PHT. However, also by using only 
the 5-fold preimage of the PHT the same MPEP was ob- 
tained. Thus neither the exact number of preimages nor 
the accurate determination of the PHT is crucial - an in- 
dication about the robustness of the method. 

In Fig. Othe MPEP is compared with a direct Monte 
Carlo simulation of the escape trajectory. The simulation 
was carried out with a noise strength of D = 2.4 x 10~^, 
resulting in a trajectory of length 2.5 x 10-'^^. One sees a 
strikingly good agreement of about 15 iterations before 
the path approaches the saddle point, reaching back to 
the 10-th preimage of the PHT. This remarkable coinci- 
dence of the simulated path with the MPEP corroborates 
the general validity of the procedure. 

To connect the MPEP with the low-period unstable 
periodic orbits embedded in the CA, we calculate all pe- 
riodic orbits up to period 20 . Only one periodic orbit 
is found to be close to the MPEP, having period 9. It 
is also included in Fig. |31 Thus, if one had adopted the 
previous method of finding periodic orbits close to the 
MPEP by using stochastic simulations of the system, one 
might have been led to conclude that the MPEP starts 
at that period 9 orbit. However, such a reasoning ob- 
scures the general deterministic structure of the MPEP 
with respect to the PHT, as presented in this work. 

Next we exemplify the method with the Ikeda map 
[iif , a discrete model of a laser pulse in an optical cavity. 
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FIG. 3: (color online). The MPEP is depicted with crosses, 
connected through a dashed line to guide the eyes. The es- 
cape path, as obtained from a Monte Carlo simulation with a 
noise strength of D = 2.4 x 10"'', is shown with squares. An 
unstable period-9 orbit is displayed with diamonds. 



With complex variables it has the form 



6z„ exp 



IK 



IT] 



1 



^n, (6) 



where z„ = Xn+iy-n is related to the amplitude and phase 
of the nth laser pulse exiting the cavity. The parameter 
a is the laser input amplitude and (1 — 6) the damping, 
while the empty cavity detuning is given by k and the 
detuning due to a nonlinear dielectric medium by rj. 

We fix the parameters at a = 0.92,6 — 0.9, k — 0.4 
and rj = 6.0. For the noisefree system, two stable states 
are present, a fixed point and a CA. In Fig. ^we present 
the result for the Ikeda map, analogously to Fig. ^ The 
CA, the basin of attraction and the saddle point on the 
boundary (square) are depicted. It is also shown the PHT 
closest to the boundary together with its 5-fold preim- 
ages (circles) and the MPEP (crosses). We include two 
period-5 orbits as well (triangles), which are located near 
the images of the PHT. The behaviour of the MPEP and 
the optimal force (not shown) is qualitatively the same 
up to where the path reaches the PHT. Then it passes 
very close to two unstable period-5 orbits, which are out- 
side the CA, before it approaches the basin boundary. 
This is so because they happen to be located exactly at 
the regions of phase space where the noise elongations 
attain a maximum, which occurs close to the images of 
the PHT 23] . The passing of the MPEP through another 
invariant set was also found in a different system [2l| , as 
well in the more complex situation when there exists a 
chaotic saddle 29]. It does not happen for the Henon 
map though, since all the unstable periodic orbits lie, for 
the present parameter values, on the CA. 

To conclude, we have unveiled the general mechanism 
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FIG. 4: (color online). CA for the Ikeda map with a — 
0.92, b — 0.9, K = 0.4, and -q — 6.0. The hatched region marks 
the basin of the fixed point located at {x,y) = (2.995, 3.947). 
Circles depict the PHT at {x, y) = (1.327, 0.559) together with 
its 5-fold preimages (numbered -1,.. ,-5). The MPEP is shown 
with crosses, connected through a dashed line to guide the 
eyes. In addition, two unstable period-5 orbits are displayed 
with upward and downward triangles, respectively. The sad- 
dle point on the boundary is marked with a square. 



of noise-induced escape from a nonhyperbolic CA. It is 
shown that the MPEP exits at the vicinity of the PHT 
closest to the basin boundary. Our findings are estab- 
lished by solving unambiguously the auxiliary Hamilto- 
nian system, which yields the exact description in the low 
noise limit. This has, for the first time, been established 
without taking recourse to Monte Carlo simulations, and 
only the knowledge of the deterministic structure of the 
PHT is required. This mechanism gives a robust prac- 
tical procedure and should also be able to be confirmed 
experimentally advantageous both for stabi- 

lizing systems and for energy-optimal switching between 
different states (cf. [ilj). 
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